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Proton structure effects in hydrogenic bound states are analyzed using nonrelativistic QED effec- 
tive field theory. Implications for the Lamb shift in muonic hydrogen are discussed. Model-dependent 
assumptions in previous analyses are isolated, and sensitivity to poorly constrained hadronic struc- 
ture in the two-photon exchange contribution is identified. 
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Introduction. Atomic spectroscopy can provide the 
most precise determination of fundamental hadron prop- 
erties, such as the proton radius [lH3|. The need for sys- 
tematic analysis to translate between bound state ener- 
gies and hadronic observables is sharpened by a discrep- 
ancy between the recent muonic hydrogen Lamb shift 
measurement Q and existing theoretical calculations. 
Using a mo del- independent extraction of the charge ra- 
dius from electron scattering data (r = 0.871(10) fm f5|; 
see also 0, Q ) or an extraction from electronic hydro- 
gen spectroscopy (r = 0.8768(69) fm Q; see also Q), the 
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gen lies 0.258(90) mcV or 0.311(63) meV above theory. 
The discrepancy brings into question the treatment of 
proton structure effects in atomic bound states, and has 
generated speculations on new forces acting in the muon- 
proton system Q , inadequate treatment of proton charge 
density correlations Q , and modifications of offshcU pho- 
ton vertices [lot . 

Non-relativistic QED (NRQED) [ill is a field theory 
describing the interactions of photons and nonrelativistic 
matter. The NRQED lagrangian is constructed to yield 
predictions valid to any fixed order in small parameters 
a and \q\/M, where |qf| denotes a typical bound state 
momentum, and M is a mass scale for the nonrelativistic 
particle. NRQED provides a rigorous framework to study 
the effects of proton structure, avoiding problems of dou- 
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blc counting in bound state energy computations [l3 |: 
eliminating difficulties of interpretation for the polariz- 
ability of a strongly interacting particle [l3| ; and provid- 
ing trivial derivations of universal properties, such as the 
low energy theorems of Compton scattering [l3 | . 

We examine the NRQED framework for determin- 
ing proton structure corrections in atomic bound states. 
The Lamb shift in muonic hydrogen is the first mea- 
surement directly sensitive to the spin-independent, pro- 
ton structure-dependent, contact interaction appearing 
in NRQED {d2 below). The strength of this interaction 
is not determined by measured on-shell form factors, or 
inelastic structure functions of the proton. We identify 
model-dependent assumptions in previous analyses and 
discuss whether poorly constrained proton structure cor- 
rections can account for the above-mentioned discrep- 
ancy. We conclude by outlining extensions of the theo- 
retical analysis and related applications. 

NRQED. Consider the formalism for electron-proton 
bound states; the substitution e — s- /i applies for the 
muon-proton system. The NRQED lagrangian can be 
decomposed as £nrqed = Cy + Ce + Cp + ^contact- 
£^ contains the photon kinetic term and vacuum po- 
larization corrections; these corrections can be treated 
separately and will not be considered here. Through 
0(1 /ml), ^M,\l^ 
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Here ipe is a two-component spinor representing the non- 
relativistic electron field, cr is the Pauli spin matrix, Dt 
and D are covariant derivatives and E, B are the elec- 
tric and magnetic fields. Prefactors are chosen for con- 
venience so that for a point-like fermion at tree level, 
cf = CD ^ cs = cwi = CAi = 1 and cw2 = Cp'p = cm = 
CA2 = 0. A similar expression holds for Cp with e — — Ze 



(Z = 1 for the proton). Relevant contact interactions in 
the single proton plus single electron sector are 

i-contact — «1 h "2 ■ (^ ) 

m^mp m^-mp 

The coefficients q, di depend on the choice of ultravi- 
olet regulator. Since no new bound state computations 
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FIG. 1: Two photon exchange amphtude for e p ^ e p scat- 
tering. 



are necessary, we will quote results for phenomenological 
inputs and bound state energies that are independent of 
this choice. We are interested in proton structure correc- 
tions to energy levels through order rrv^a^ / m^, and there- 
fore need cf,d,s in through C(a), and di,2 through 
0(0^). Other operators in Cp will enter when we analyze 
the low expansion of the forward Compton amplitude 
to constrain d2- Knowledge of the c^'s and d^'s allows us 
to determine corrections to energy levels. For example, 
coefficients 0^°*"°"^ and c?2 lead to first order energy shifts 
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where rrij. = 'm^mp/{m^ + rup) is the reduced mass. 
Matching. The NRQED Wilson coefficients are de- 
termined by enforcing matching conditions between full 
and effective theories using convenient low energy observ- 
ables. We concentrate on the matching conditions for the 
proton. 

One photon matching. Wilson coefficients for operators 
coupling to a single photon are determined in terms of the 
proton elastic form factors and their derivatives at = 
by using ([T|) to compute the amplitude for elastic scatter- 
ing of a proton via the electromagnetic current [IB-Ill- 
The form factors satisfy Fi(0) — 1, F2{0) = Bp, 
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where A is a photon mass [l9[. These expressions serve 
to define the phenomenological parameters ap « 1.793, 



and r^t 
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Two photon matching. The coefficients cai, cai can be 
determined by comparing to spin-averaged amplitudes 
for forward and backward Compton scattering in the lab 
frame [17 1 
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12.0(6) X 10-''fm^ and /3 = 1.9(5) 



where a = 

10- W [201 . 

Contact interactions. The coefficients in Q can be fixed 
using the zero- momentum limit for e~p^e~p scattering, 
cf. Fig.[TJ The tree level, 0{a), amplitude is reproduced 
by the effective field theory, and the d^'s receive a non- 
zero contribution starting at 0{a^). We focus on the 
spin-independent case and neglect higher order radiative 
corrections. The relevant proton matrix element is the 
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Our normalizations are such that for a point parti- 
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The matching condition for the spin-averaged zero- 
momentum four-point amplitude is 
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where 5d2 denotes the contribution to ^2 in addition to 
the point particle value. 

The imaginary part of the WiS can be related to mea- 
sured quantities. By inserting a complete set of states 
into ^ , the proton contribution to Im Wi is expressed 
in terms of proton form factors, and the continuum con- 



tribution to ImWi is determined by inelastic structure 
functions. Using dispersion relations, W2 can be fully re- 
constructed from its imaginary part. Since Wi requires a 
subtraction for a convergent dispersion relation, knowl- 
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edge of Wi (0, Q^) is also needed. Thus Wi can be written 

W2{v, Q') = M/f "(i., + W^^^iy, , (8) 

where the superscript numbers denote the number of sub- 
tractions. The proton terms are 



with Fi = Fi{—Q^). The continuum terms are 
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where i'cut(g ) = g + imT^mp + is the threshold 
for pion production, and ImWi(i', g^) are proportional 
to inelastic scattering cross sections. 

The small and large limits of Wi(0, g^) can be 
studied in a model-independent way. Using NRQED to 
compute the amplitude for double scatterin g of a proton 
in an external static magnetic field we find hj 

W^i(0,g2) ^ 2ap{2 + ap)+Q^{2mpP/a-ap/ml 
(2/3) [(1 + apfirl^r - (r^)'] } + OiQ^) . (11) 

At g^ ;» GeV^ we may evaluate ^ using the opera- 
tor product expansion (OPE). Leading terms arise from 
dimension-four operators and scale as [l3|. The in- 
termediate g^ region is not constrained by existing mea- 
surements. This lack of knowledge about Wi(0,g^) in- 
troduces model-dependence in the theoretical prediction 
for the Lamb shift, which has so far been ignored in the 
literature. 

Given this model-dependence, how did previous stud- 
ies obtain numerical predictions? The most common ap- 
proach is to pretend that 1^1(0, g^) can be separated 
into "proton" and "non-proton" contributions. A "pro- 
ton" part for Wi(0, Q^) is obtained by inserting the ver- 
tex with onshell form factors into the Feynman diagrams 
for a relativistic pointlike particle. For definiteness we 
refer to this approach as the "Sticking In Form Factors" 
(SIFF) model. Explicitly, 
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(0,g2) = 2F2(2Fi+F2), 
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We emphasize that ([T^ is not derived from a well defined 
local field theory. In fact, no local lagrangian can give 
such Feynman rules. Note also that Wf^^^(0,g^) does 
not have the correct large behavior. A "non-proton" 
part is obtained by multiplying the 2mpQ^P/a term in 



by a function of Q [21|, |23 . The models used in 



[2ll |22| again do not have the correct large behavior. 



Unlike Im Wi, we stress that the separation of Wi(0, g^) 
into proton and non-proton parts is not well-defined. 
Bound state energies. The use of an effective field 
theory allows us to systematically classify the proton 
structure corrections to energy levels. Using ([3]), proton 
vertex corrections, of order {Za)^ and {Z'^a){Za)*, are 
determined by cu. Our definition (j4|) of the proton ra- 
dius in the presence of radiative corrections implies that 
the Z'^a{Za)'^ correction is unchanged from the point- 
particle result, so that 17 
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Two-photon exchange corrections, of order (Za)^, are 
determined by c?2. Considering (|8]), it is natural to de- 
compose the correction as 



~^ = (5£:P''°t°n -|- ^£;continuum _|_ (0,Q^) ^ (•^4) 



It is convenient to subtract limQ2^Q W^i(0, g^) from 
1^1(0, g2) in ©, and add \im.Q2^oWi{0,Q^) to 
iyf^(i/,g2). Infrared singular terms in ([7|) are then 
confined to the proton pole contribution. Having fixed 
this terminology, we proceed to discuss each of the three 
terms in SE^'"°^'^ in turn. Our discussion so far applies to 
general hydrogenic bound states. To investigate numeri- 
cal results we now specialize to muonic hydrogen ( "/i_ff " ) . 
Proton pole contribution. We content ourselves with a 
simple dipole model for the elastic form factors. 



GM{q')/GM{0)^[l-qVAT\ (15) 



where Ge = Fi + {q^/4m^)F2, Cm = Fi + F2 and 
A^ = 0.71 GeV^. Wc return to a more sophisticated anal- 
ysis of this contribution, and analo gou s spin-dependent 
contributions, in forthcoming work 17[. After isolating 
the finite term in ([7]), for muonic hydrogen 



,5£;P^*°"(n5) « (8/n3) (0.016 meV) 
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We refrain from giving a detailed error estimation here; 
for the purpose of explaining the muonic hydrogen 
anomaly, an error smaller than 100 % does not have very 
substantial impact. 

Continuum contribution. A recent determination of the 
continuum contribution is 12111 



se: 



continuum 



(nS) « (S/n^) (-0.0127(5) meV) , (17) 
0.014(2) mcV [12, 



in line with pre vious results, 
-0.016(3) meV [23. 
Wi{0,Q^) contribution. In the SIFF model ^ one 
finds, 



SE 
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{nS) = (S/n^) (-0.034meV) . (18) 
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n c fool 

Her. 1221 


Ref. [25] 


This work 


vertex 


-0.0099 


-0.0096 


-0.0108 
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0.035 


0.051 


Model Dependent 




r j-'eontinuum 






0.013 [21] 


Total 


0.025 


0.042 





TABLE I: Comparison between this and previous works for 
O(a^) proton structure corrections to the 2P — 2S Lamb shift 
in muonic hydrogen, in meV. 



The sum of the proton pole and W^i(0, Q"^) contributions 
in this model, 0.016 meV - 0.034 meV = -0.018 meV, 
reproduces previous results [131 • It is not hard to con- 
struct model functions for Wi{0,Q^) that have the cor- 
rect small-Q^ and large-Q^ behavior, but give a much 
larger contribution than the SIFF model. 
Comparison to previous results. In order to make the 
comparison to the literature clearer, we collect the re- 
sults of this analysis in Table HI which compares numeri- 
cal results for 0{a^) proton structure corrections in the 
2P — 2S Lamb shift of muonic hydrogen. We focus on 
the two reference sources that were used in namely 

[H and 0. These works model (JS^^"''^'^ as a sum 
of proton and non-proton contributions, adding the re- 
spective terms to SE^'°^°'' and J^j^^tinuum^ ^^^^^ 
simplify the comparison we present in the table the total 
contribution to 5E^^^~'^ from [g^j and [2^. In partic- 
ular for [2^ we add the {Za)^ nuclear size correction 
(0.0232 meV) and the proton polarizability correction 
(0.012 meV). For ^ we add the {Zaf nuclear size cor- 
rection (0.0232 meV), the polarizability correction (0.015 
meV), and the recoil finite size correction (0.013 meV). 
In [21 th e nuclear size correction at order {ZaY' from [22| 
andjii employs the SIFF ansatz ^ for VFi(0, Q^); the 



(r^)'^ scaling employed in Q assumes the large rrip limit 
and a one-parameter model for Ge and Gm- 

Let us note three differences between our results and 
the theoretical predictions used in Q, and collected in 
Table [B First, the proton vertex correction from 
22I I24II uses a different convention for the c harg e ra- 
dius [26[ , while the result from 27 1 , adopted in 25 1 , uses 
a model-dependent SIFF prescription for the proton ver- 
tex correction; the complete result with the charge ra- 
dius definition (j4]) is given by ([T3|) . displayed in the first 
line of the table. Second, the "recoil finite size" of [25| . 
adopted from is in fact part oi5E^^^~'^ . Including it 
as separate contribution would lead to double counting. 
Third, the SE^^'^^''^ contribution is model-dependent; 
the current theoretical prediction is based on the SIFF 
ansatz. We conclude that the dominant radiative correc- 
tion to proton structure is subject to uncertainties from 
unreliable hadronic models. 

Discussion. We have presented the NRQED formalism 



for systematically analyzing proton structure effects in 
hydrogenic bound states. The Lamb shift in muonic hy- 
drogen is sensitive to a new structure-dependent contact 
interaction The strength of this interaction is not 
determined by measured proton form factors or inelas- 
tic structure functions. Taking all other contributions 
as fixed, the muonic hydrogen Lamb shift determines d2 
in ([2]). NRQED then predicts a universal shift for other 
spin-independent energy splittings in muonic hydrogen. 

The strength of the contact interaction can be related 
to a so-far poorly constrained piece of the forward Comp- 
ton amplitude of the proton, Wi{Q,Q^). In this Letter, 
we have established some mo del- independent properties 
of W^i(0,Q2). Firstly, the O^Q"^) Taylor expansion ^ 
is determined by NRQED in terms of measured quanti- 
ties; secondly, the asymptotic behavior is determined by 
OPE techniques to be ~ . The intermediate region 
remains poorly constrained [2^. The lack of theoretical 
control over Wi{Q,Q^) introduces theoretical uncertain- 
ties that have not been taken into account in the liter- 
ature. A common approach is to use the SIFF model 
([T^. but this is not derived from first principles and 
gives the misleading impression that the dominant de- 
pendence is constrained by onshell form factors [U, [g^l ■ 
Such extrapolations represent models for VFi(0, Q^), typ- 
ically without the correct large behavior. While we 
do not attempt an exphcit modeling of Wi(0, Q^), we 
believe that the uncertainty assigned to this contribu- 
tion (< 0.004 me V 0) is underestimated by at least an 
order of magnitude. 

As further applications, the nonrelativistic effective 
theory for vector fields can similarly be employed to 
describe deuterium. The NRQED lagrangian at order 
1/M^ can be used to systematically analyze (5i?*'"°~T in 
the small-lepton mass limit relevant to electronic hydro- 
gen, and describes spin polarizabilities of the proton [17 1 . 
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